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Solving equations (3) and (4) for y(s) and 0(s) gives the Laplace 
transform of y(x) and 0(x), respectively, as 



where 

q(s) = s 6 + 6s^ - os 2 + 70 - oB 

Goland and Luke (reference 4) showed that y(x) and 0(x) could he 
written as a converging series by expanding the transforms (5) and (6) 
in terms of symmetric polynomials of the squares of the roots of q(s) 
and applying the inverse transform. A discussion of this expansion is 
given in section 4 of appendix A where it is shown that l/q(s) can be 
written as 



where 

M 



Tg = 5^ + a 

T = -63 - as -p 7 
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For n J> 3 , 


T n = -&V! + al n _2 + (aB - P7)T n _3 


(8) 


When the series expansion of l/q(s), equation (7) I, Is substituted 
. + on-nations (b) and (6) the transforms y(s) and 0(s) Become sums of 
infinite series with terms of two distinct types J that Is, terms of types 


A_ 

s* 


and 

Be~ SX ° 

gm 

where m is a positive integer* 

The inverse Laplace transform of ~ (see pair no. 3, P‘ 2 95, of 
reference 9 ) for x > 0 is 



(9) 


and the Inverse Laplace transform of 
of reference 9) for x >Xq^ 0 is 



s m 


(see pair no. 63, P* 298 , 



B(x - xo) 

O - T) 


m- 

T" 
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FLUTTER OF A UNIFORM.. WING WITH AN ARBITRARILY 
PLACED MASS ACCORDING TO A DIFFERENTIAL- 
EQUATION ANALYSIS AND A COMPARISON 
WITH EXPERIMENT 

By Harry L. Runyan and Charles E. Watkins 


SUMMARY 


A method is presented for the calculation of the flutter speed of a 
uniform wing carrying an arbitrarily placed concentrated mass. The 
method, an extension of recently published work by Goland and Luke, 
involves the solution of the differential equations of motion of the 
wing at flutter speed and therefore does not require the assumption of 
specific normal modes of vibration. The order of the flutter determi- 
nant to be solved by this method depends upon the order of the system 
of differential equations and not upon the number of modes of vibration 
involved. 

The differential equations are solved by operational methods and 
a brief discussion of operational methods as applied to boundary- value 
problems is included in one of two appendixes . A comparison is made 
with experiment for a wing with a large eccentrically mounted weight 
and good agreement is obtained. Sample calculations are presented to 
illustrate the method j and curves of amplitudes of displacement, torque, 
and shear for a particular case are compared with corresponding curves 
computed from the first uncoupled normal modes. 

For convenience, the method employs two-dimensional air forces 
and could be extended to apply to uniform wings with any number of 
arbitrarily placed concentrated weights, one of which might be considered 
as a fuselage. The location of such masses as engines, fuel tanks, and 
landing-gear installations might be used to advantage in increasing the 
flutter speed of a given wing. 


INTRODUCTION 


The common procedures in flutter analysis of an airplane wing 
involve many simplifying assumptions. In particular the degrees of 
freedom of the wing are usually determined by choosing the first few 
normal modes of the structure, and the wing motion at flutter is then 
described in terms of these chosen modes. This approach of employing 
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prescribed inodes is often adapted to the Rayleigh type analysis of 
vibration and may he referred to as Rayleigh type analysis . In 
specific calculations with this method the amount of work required is 
proportional to the number of normal modes involved. In particular, the 
order of the flutter determinant that must be solved depends directly 
upon the number of modes involved. For simple wings, without concen- 
trated masses, the Rayleigh type analysis usually yieldB satisfactory 
results with not more than two or three normal modes. However, if the 
wing carries concentrated masses, such as engine, fuel tank, or landing- 
gear installations, so many normal modes may be required to obtain satis- 
factory results that the Rayleigh method may not be the most feasible 
method . 

Tn cases where many degrees of freedom are involved the most logical 
procedure would be to treat the system of differential equations of 
motion of the wing rather than to choose specific modes. This method is 
in general very difficult and tedious to oarry through, although it has 
the advantage that the order of the flutter determinant that must be 
solved depends only upon the order of the system of differential equations 
and not upon the number of modes of vibration involved. 

As early as 1 929 Kussner (reference 1) used the differential 
equation approach to formulate the problem in the form of an integro- 
differential equation for a wing of general plan form . Kussner set up 
some particular examples and suggested a method of solution by a process 
of iteration. This method was not followed up until during the war 
when some related work was undertaken in Germany but not finished. 

Wielandt (reference 2) has recently made contributions to the treatment 
of nonself-ad joint differential equations by iterative processes. In 
the light of these contributions perhaps the problem of flutter analysis 
as proposed by Kussner warrants further investigation. 

Recently, Goland (reference 3) applied the differential- equation 
method to a uniform cantilever wing and was able to carry out the 
solution of the flutter problem by straightforward methods. In refer- 
ence 4 Goland and Luke extended the solution of the problem of the 
uniform wing to include a uniform wing carrying a fuselage at the 
semispan and concentrated weights at the tips. Goland and Luke made 
use of the Laplace transform to solve the differential equations by 
operational methods for both the symmetric and antisymmetric types of 
flutter. In both references 3 ani 4, the objective was to compare 
flutter speeds and certain flutter parameters for specific uniform 
wings calculated by the differential- equations method with the same 
quantities calculated by the Rayleigh method when only the fundamental 
bending and torsion modes were used in the calculations. Fairly close 
agreement between results calculated by the two methods were obtained 
in both references 3 and 4. No comparison with experiment, however, 
was made in either case. 

The results of a systematic series of flutter tests made to 
determine the effect of concentrated weights and concentrated weight 
positions on the flutter speed of a uniform cantilever wing are reported 
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in reference 5* After these experiments were finished, an attempt was made 
to compare the results with a theoretical analysis ty the Rayleigh method. 
In cases where the mass of the weight was of the same order as that of the 
wing and placed so that the distance "between its center of gravity and the 
elastic axis of the wing was a considerable fraction of the wing chord, 
however, several normal modes would have to "be employed and there was no 
way of knowing in advance just what number should "be used. Because of 
this difficulty and "because the wing was a uniform wing, the most extreme 
case was chosen from reference 5 and investigated "by the differential- 
equations method by following an extended procedure of Goland and Luke. 

The purpose of this paper is to report the results of this investigation. 

The paper consists of the main text and two appendixes . In the 
main text the differential-equation method is set up for any uniform 
cantilever wing with an arbitrarily placed concentrated weight and the 
solution, based on an extension of the method used by Goland and Luke, is 
developed. Application Is then made to a particular wing- weight system 
used in reference 5* and comparison with experimental results is given. 

The mass of the weight (weight labeled Ja in reference 5) was about 
92 percent of the mass of the wing and at each spanwise weight position 
the weight was placed so that its center of gravity was about 0.4l chord 
forward of the elastic axis of the wing, (it may be mentioned for the 
sake of comparison that in the numerical example treated in reference 4, 
the mass of the weight was only 39 percent of the mass of the wing and 
placed 0.1 chord behind the elastic axis of the wing.) The geometric 
aspect ratio of the wing was 6, which was considered large enough to 
warrant the use of two-dimensional air forces without aspect-ratio 
corrections for oscillating instability (not necessarily so for the 
divergent type of instability (see reference 6)). One other simpli- 
fication was the omission of terms due to structural damping. The 
computed results agree remarkably well with experimental results, 
particularly in regard to trends. 

In appendix A the method used by Goland and Luke, which includes 
the derivation of the differential equations, for a wing carrying a tip 
weight is outlined and extended to a wing carrying an arbitrarily placed 
weight. A somewhat general but brief discussion of operational methods 
of solving boundary-value problems is included and illustrated with a 
simple example for readers who might be interested but are not familiar 
with the operational approach. 

In appendix B the derivation of the flutter determinant Is com- 
pleted and a method of solving the determinant Is illustrated by a 
detailed calculation of the flutter speed for the wing and one weight 
position of the wing-weight combination discussed in the test. As a 
final topic in this appendix the solution obtained for the flutter 
determinant is used with the solutions of the differential equations to 
calculate the amplitudes and phase angles of the deflection curves of 
the wing- weight system at flutter speed. 
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simbois 


(The symbols are given in terms of a consistent set of units that 
are convenient for the computations in this paper. They can he converted 
to any desired set of units by proper attention to the dimensions 
involved.) 


nondimens ional distance of elastic axis from midchord 
measured in half- chords, positive for positions of 
elastic axis behind midchord 


b wing half- chord, feet 

©1 chordwise distance of wing center of gravity from 

elastic axis, positive for center of gravity behind 
elastic axis, feet 

©2 chordwise distance of weight center of gravity from 

elastic axis, positive for center of gravity behind 
elastic axis, feet 


8 

I 


w 


gravitational constant, feet per second per second 

mass moment of inertia of uniform wing per unit of 
spanwise length, referred to wing elastic axis. 


pound- second 


MK- 


!) 


k 2 

k 

L 

h + 

l 

Z 1 


mas3 moment of inertia of weight referred to wing 

p 

elastic axis, foot- pound- second 

radius of gyration of wing sections about wing elastic 
axis, feet 

radius of gyration of weight about elastic axis, feet 

reduced- frequency parameter ( ~ J 

aerodynamic lift force per unit of spanwise length 


p 

jtpb Lh 


w» 3 [L a - *1,(2 + a )j 

semispan of ving, feet 

location of weight measured from wing root, feet 
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aerodynamic coefficients as tabulated in reference 7 
M 


aerodynamic moment per unit of spanwise length taken 
about elastic axis 


M- 


[y + iMy' = Jtpb^ Jmji - L h (| + a)] 

Mq + IMq 1 = TtphA 


m 

*w 

N 

T 

e 1 ^ 2^3 

s 

t 

T 

x n 


y 

tcs 

x 

Y(x,t) 

y(*) 

EIo. 


M a - M 2 + “ 


- M h| + ») + Lh(| + «J 


weight of wing model , pounds 

mass of wing per unit length 

weight of concentrated weight , pounds 

transverse shear force in wing at station x 

torsional moment in wing at station x 

roots of cubic equation 

operator used in Laplace transformation 

time coordinate 

sum of all symmetric polynomial functions in K^_, Tfej ^3 
which are of degree n 

experimental flutter speed for wing without weight, 
feet per second 

flutter speed, feet per second 
reduced flutter speed 

spanwise coordinate measured from wing root 

general mode shape function in tending 

mode shape function in tending after assumption of 
harmonic motion ^y^(x) + iy^(x)^ 

2 

flexural rigidity of uniform wing, pound-feet 


GJ 


torsional rigidity of uniform wing, pound- feet^ 
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a - C( B + ^ + 

+ Lg + ILg y 



5 “ sK 1 + “e + “V 


p 

A 

0(x,t) 

0(x) 

03 


mass ratio 


rtpV; 

m 


air density, slugs per cubic foot 
complex value of determinant 

value of A when real and imaginary parts are equal 

general mode shape function in torsion 

mode shape function in torsion after assumption of 
harmonic motion (©2(x) + i 03 (x)j 

circular frequency at flutter, radians per second 

frequency, cycles per Becond 


ANALYSIS 


As mentioned in the introduction the differential equations that 
govern the motion of a uniform wing at flutter speed, as derived by 
Goland in reference 3, and a method of solving the equations for a 
uniform cantilever wing carrying an arbitrarily placed weight, based on 
a method developed by Goland and Luke in reference 4, are discussed in 
appendix A. This section, therefore, is devoted to a brief discussion 
of the differential equations of motion of the wing, the boundary condi- 
tions, solution of the boundary-value problem by means of the Laplace 
transform, and the solution of the flutter determinant. 

The differential equations and boundary conditions that govern 
the motion, at flutter speed, of a cantilever wing of length 2 with 
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a concentrated weight placed Zp units along the span from the root 
section and e2 units forward of the elastic axis of the wing, as derived 
in appendix A, are 


y lv (x) - ay(x) - (30 (x) = 0 (l) 

0 ' T (oc) + 7 y(x) + 50(x) = 0 (2) 


(a) y(o) = y'(o) = 0 ( 0 ) = 0 

(h) Elpy'^Z) =EI b y'"(z) = GJ0 ' ( Z) =0 

(c) EI- b [y"’(Z 1 - 0) - y ,M (Z x + 0)J = " + e 2 0 ( 2 l)] 

(d) . GjjV (Zp - 0 ) - 0 T (Zp + 0)>^[e 2 y(Zi) + K 2 2 0 (ip)] 

where 

a = f^( m + b + 1 V) 

P = + Lq + iLe y 

T " + “y + “V') 

6 = s£(i + Me + ate y 


and where y(x) is the displacement of a chordwise element of the elastic 
axis of the wing at span position x due to "bending j 0(x) is the corre- 
sponding displacement due to torsionj primes associated with y and 9 
Indicate differentiation with respect to xj EI b is the ^flexural rigidity 


W 


w 


of the wingj GJ is the torsional rigidity of the wingj 

weightj m is mass per unit length of wingj and to is the circular 
frequency of tending and torsion at flutter. 


is mass of the 


notation y M 1 ( Zp 
it approaches as 
and y" '(Zp + 0) 
approaches as x - 


In condition (c) the 

0) indicates that y M, (x) is to have the value that 
x — >Zp from the inboard side of the weight 


indicates that y' M (x) is to have the value that it 
Z-i from the outboard side of the weight. Similar 

_ v _ - 1 / _ « \ 


meanings are given to $*(^1 ~ 0) and 0 ? (2i + 0) 
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The quantities Ly + iLy 1 , Lg + iLg 1 , My + iMy 1 } and Mg + iMg 1 
can be written In terms of tabulated quantities as follows: 

Ly + lLy = Itpb Ljj 

Le + iLg' = npb 3 |L a - LhQ; + ajj 

My + iMy 1 = rtpb^j^Mh - Lh (7; + a)^j 

^ + «%' = "etpa - 41 + a )- 41 + ») * Jill + _ 

In reference 7 the values of L^, L a , and are expressed in terms 

of Theodorsen*s F and G functions of reference 8 and tabulated for 
various values of the reduced speed v/bco. 

The root conditions (a) and the boundary conditions (b), of the 
boundary- value problem, are the usual conditions that must be imposed 
upon the equations of a vibrating cantilever beam (or wing). Condi- 
tions (c) and (d) stipulate discontinuities of determinable magnitudes 
in transverse shearing force and torque, respectively* 

Applying the Laplace transfoim (see appendix A) 

0 

e~ sx f(x) cbt = f (s) 

to equations (l) and (2) and making use of conditions (a), (c), and (d) 
gives 

- sY 2 - Y3 + e _s2 lJ7"'(Zi - 0) - y'"(Zi + 0)J - ay(s) - ( 30 (s) = 0 

( 3 ) 

and 

s^6(s) - 0i + e S 1 ( 2i - 0) - d'(l± + 0)J + 80 (s) + 7y(s) = 0 ( 4 ) 

where 

y 2 = y "(o) 

Y 3 = y ,M (o) 

0 1 = e , (o) 
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and 


-2L m 2n+l 


2n+4 


2n+5 


V“ T n x V — T n x T n x 

•W ■ e ^2_ (^TljT ' ^LiSTTmT - 7I 3 2_(5^)T 

n=0 n=0 n=0 


- Spa 


r^a 


J&1+5 


i—- (2n + 5) ! 
n=0 


e'Ui - 0) - e'Ui + 0) 


n=0 


TnU - 

(2n + 5) ' 


4 ® ^ + D : - 



y'" 




(12) 


where In both equation (ll) and equation (12) the terms involving (x 
are to "be considered as zero when x = l]_. 


*l) 


Equations (ll) and (l2) are general expressions for the amplitudes 
or displacements of a point x of the elastic axiB of a uniform wing 
vibrating in bending and torsion under the conditions of flutter with an 
arbitrarily placed concentrated weight. When the weight is concentrated 
at the wing tip the equations correspond to those obtained by Goland 
except for a difference in root conditions. When the weight is con- 
centrated at the root (or if the mass of the weight is reduced to zero) 
the equations reduce to those for a uniform cantilever wing. These 
equations may appear rather formidable in their present formj however, 
only the first few terms of each summation seem necessary for most cases. 


In the derivation of the flutter determinant in appendix B it is 
shown that since terms involving (x - Zp) drop out of both equation (ll) 
and equation (12) at x = 2p, the values of y(Z^) and 0(2^) can 
obtained from the terms not involving (x - 2]). Then, by making use of 
conditions (c) and (d) again, linear expression in Y2, Y3, and 6^^ can 
be substituted for the bracketed expressions 


- 0) - 
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For n J> 3, 


T n = " &T n-l + <^-2 + (aB - Pr)T n -3 (8) 

When the series expansion of l/q(s), equation (j ) , is substituted 
into equations ( 5 ) and (6) the transforms y(s) and 0 (b) become sums of 
infinite series with terms of two distinct types } that is, terms of types 

A_ 

s* 1 

and 

Be~ sx ° 

gm 

where m is a positive integer. 

The inverse Laplace transform of — (see pair no. 3, p. 295, of 
reference 9 ) for x > 0 is 



Ax™' 1 


(9) 


B ~ BX o 

and the inverse Laplace transform of 6 _ — (see pair no. 63 , p. 298 , 


of reference 9) for x > Xq 0 is 


gm 


r-lJSe, 


3I o 


0m 


B(x - xo)”' 1 
(m - l) ! 




( 10 ) 
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When the egression for l/q(s) from equation (7) is substituted into 
equations (5) and (6) and the inverse transforms is applied, the following 
series expressions of y(x) and 6(x) can be obtained: 


y(x) = Y c 


00 

2n+2 00 

2n+4 


00 2n+5 

y - 

T n* , ,V 

T n x 

+ Y 

5 V T n x 

Z_ 

(2n + 2) J + 5 / 

(2n + 4) I 

*3 

V (2n + 5 ) i 

__ n=0 

n=0' 

— 


- n=0 


00 _ m 2n+3 
x n JL 


n=0 


(2n + 3) I 


V- T n x 2n+ 5 

1 Z_ (2n + 5) l 
n=0 


- P 


e'Ui - 0) - e'(ii + 0) 


T n (x - Z l} 2n ^ 
(2n + 5) ! 


y ,M (ll - 0) - y>"(2i + 0) 


n=0 


T n (x - Zp) 211 ^ 
(2n + 5) l 


n=0 


T n (* - l l> 

(2n + 3) I 


2n+3 


(n) 
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and 


V 24 T n x 2n+1 T n / n+4 

•w - e l > , * ~ - ^2 > — — — - rc 


(2n + 1) ! 

n=0 


_!L n, 2n+5 


(2n + 4)1 3 Z (2n + 5)! 

n=0 n=0 


- 0ia 


T n x 2n+ 5 


V 

4 — ( 2 n + 5 ) i 
n=0 


o'(Zi - 0) - e'Oi + 0) 


a 


-l L n=0 


rnU - 1^2 

(2n + 5) I 


r n( x - Z l) 2n+1 


n=0 


(2n + 1) I 


+ 7 


y '"0i ■ °) " y '"0i + °) 


y T n(* - hf^Z 
4s) ( 2n + 5 ) 1 


( 12 ) 


where in "both equation (ll) and equation (12) the terms involving (x - ZjJ 
are to he considered as zero when x = 2]_. 


Equations (ll) and (l2) are general expressions for the amplitudes 
or displacement of a point x of the elastic axis of a uniform wing 
vibrating in bending and torsion under the conditions of flutter with an 
arbitrarily placed concentrated weight. When the weight is concentrated 
at the wing tip the equations correspond to those obtained by Goland 
except for a difference in root conditions. When the weight is con- 
centrated at the root (or if the mass of the weight is reduced to zero) 
the equations reduce to those for a uniform, cantilever wing. These 
equations may appear rather formidable in their present formj however, 
only the first few terms of each summation seem necessary for most cases. 


In the derivation of the flutter determinant in appendix B it is 
shown that since terms involving (x - Zp) drop out of both equation (ll) 
and equation (l2) at x = 2^, the values of y(Z^) and 0(2^) can be 
obtained from the terms not involving (x - 2 2). Then, by making use of 
conditions (c) and (d) again, linear expression in Yp, Y3, and 0-^ can 
be substituted for the bracketed expressions 


- 0) - y'"(u + 0) 
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and 


”e’(*i - °) - 6 '( ! i + °)] 

After the substitutions are made, equations (ll) and (12) will contain 
only the three undetermined coefficients Yg, Y 3 , and 0 -j_ for any 

particular wing-weight system of the type under consideration. Observe 
that conditions (b) have not yet been used. If these conditions are 
now imposed upon the equations, there is obtained a system of three 
linear homogeneous equations in Y 2 , Y 3 , and 0p that may be written 
for reference as 


AjY2 + BiY 3 + Cj0i = 0 


(13) 


where i = 1 , 2 , and 3 - 

The condition that a system of equations such as equations ( 13 ) have 
solutions other than the trivial solution 


Y 2 = Y 3 = 0i = 0 


Is that the determinant of the coefficients A± } B±, and Ci vanish 
(reference 10) . This corresponds to the border-line condition between 
damped (stable) and undamped (unstable) oscillations or to the point at 
which flutter occurs. It will be noted that the order of this determi- 
nant depends only on the order of the system of differential equations. 

The actual coefficients corresponding to A^, B^, and are 

complex functions of the frequency u>, the reduced flutter speed v/bco, 
and certain determinable characteristics of the wing-weight system. The 
true flutter speed is easily calculated when corresponding values of cd 
and v/bx are known. These quantities may therefore be considered as 
(the only) variable parameters in the determinant of coefficients and 
the problem of finding the true flutter speed Is reduced to that of 
finding corresponding values of these parameters that cause the determi- 
nant, hereinafter called the flutter determinant, to vanish. If v is 
set equal to zero the air forces drop out and the resulting determinant 
gives the coupled modes of vibration of the wing in still air. On the 
other hand, if u> is set equal to zero the nonoscillatory or divergence 
condition is obtained. 

Several ways of solving the flutter determinant are mentioned in 
reference 6 . Although more informative methods exist, a graphical method 
was adopted for the present work. For example, a value is assigned to 
one parameter, preferably v/txsj the flutter deteiminant is then evaluated 
for this value of v/bo and several valueB of the other parameter cd. 

The values of the flutter determinant obtained in this manner are complex 
numbers and if the real and imaginary parts of a sufficient number of 
determinant values are separauely plotted against n, tie point or points 
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■where the real and Ima ginary parte are equal are obtained. If this 
process for other values of v/bo> is repeated, a locus of determinant 
values with equal real and imaginary parts can he plotted against 
both v/bcu and 00. When enough points are determined these plots give 
the values of v/bto and cu that cause the determinant to vanish. 

An illustration of the process of solving the flutter determinant 
as described in the preceding paragraph is given in appendix B which 
contains the complete solution of the determinant for one weight position 
of the particular wing-weight system desci’ibed in the section entitled 
Application to a Specific Wing-Weight System. In general, when solving 
the flutter determinant by the preceding method, if the assumed values 
of v/bco and <o are in the neighborhood of their true values, only a few 
points need be computed to obtain a solution. In the absence of experimental 
values of these parameters and In view of the work involved In determining 
other parameters that depend on v/bto, it will be found advisable to use 
simplified methods to obtain approximate values with which to start the 
solution. 


APPLICATION TO A SPECIFIC WING-WEIGHT SYSTEM 


Attention is now turned to the application of the boundary- value 
problem discussed in the foregoing section to a specific problem. The 
wing-weight system that has been analyzed consists of a particular 
unifonrt cantilever wing and weight combination described in reference 5* 
The weight was considered as concentrated at different specified span 
positions but always at about 0 .41 chord forward of the elastic axis 
of the wing. This weight was selected because of its high' mass compared 
to that of the wing and because of the large eccentricity due to the 
distance between its center of gravity and the elastic axis of the wing. 
Furthermore, by using only the fundamental modes, first bending and first 
torsion, the Rayleigh type analysis had failed to give any reasonable 
results for this particular wing-weight combination. Pertinent data, 
based on measured characteristics of the wing as taken from reference 5, 
with the units in feet and pounds are 


Chord, feet 2/3 

Length,, feet ........... 4 

Aspect ratio (geometric) .. - 6 

Taper ratio .......... 1 

Airfoil section NACA 16-010 

W, pounds 3*48 

I, pound- second^ 0 .00080 

ETb, pound- feet^ 977 «08 

GJ, pound-feet^ 480.56 

1/k (standard air, no weight) 32*6 

e^, feet . 0.013 
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and, "based on measured characteristics of the weight, are 


~W V , pounds . . . » 3 .182 

e2, feet -0.2728 

I w , foot- pound- second^ O.OI3625 


Calculation of the flutter parameters have "been made for the wing 
without the weight and for the wing with the weight at seven different 
positions. The calculated results are compared with experimental results 
in figure 1 and in the following table: 




Calc ulated 


Experimental 

l l 

(in.) 

e 2 

(ft) 

f 

(cps) 

v/bcjo 

V' 

(fps) 

f 

(cps) 

v/bao 

V 

(fps) 

0 


25.27 

6.29 

333 

22.1 

7.22 

334 

11 

-0.2728 

19.23 

8.23 

331 

17.4 

8.88 

324 

17 

- .2728 

28.04 

6.93 

407 

^ 6 .8 

6.81 

382 

30 

- .2728 

30.68 

8.18 

526 

(b) 

(h) 

(h) 

45 

- .2728 

25.67 

7-45 

401 

(t) 

(h) 

(b) 

46 

- .2728 

24.87 

7 .06 

368 

21.8 

8 .06 

368 

48 

- .2728 

23.60 

6.07 

300 

21.4 

7.14 

320 


a It is found in reference 5 that good flutter records for this wing-weight 
system were obtained for several spanwise weight positions between the 
root section and a point 17 inches from the root section} but with the 
weight at 17 inches from the root section the wing appeared to diverge. 
However, the oscillograph records for this case showed two possible flutter 
points, one corresponding to a frequency of l 6«3 cps and another corre- 
sponding to a frequency of 26.8 cps (only the first of these is recorded 
in reference 5) * When the weight was moved farther outward from this 
point, definite divergence was noted until the weight was at a point 
46 inches from the root section. At this point and from this point to 
the tip good flutter records were obtained. 

^Divergence . 

It will be noted in the table that all the calculated flutter 
speeds are within 7 percent of the experimental values and the calculated 
frequencies and reduced speeds are within 15 percent of the experimental 
values. The calculated flutter speeds are generally slightly higher than 
the experimental values for 2 ^ ^ 17 and slightly lower for 2 ^ > 46 . 

There is no such consistent trend in the other parameters. 

In figure 1 the ratio of both calculated and experimental flutter 
speeds for the wing with a weight to the flutter speed of the wing with- 
out a weight are plotted against span position of the weight. The 
important thing to note in examining figure 1 is that the shape of the 
theoretical curve follows the shape of the experimental curve very 
closely in the regions where experimental flutter was obtained. The 
horizontal dashed line in figure 1 represents the divergence speed for 
the wing as computed by the method of reference 11 . Although the 
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correct divergence speed for different weight positions would probably 
vary, being somewhat lower with the weight at the tip than at the root, 
owing to the effect of the presence of the weight on aerodynamic forces, 
the agreement of the approximate value with experimental values is 
satisfactory . 

General expressions for the deflection curves are derived in 
appendix B from which amplitudes and phase angles for curves of deflection, 
slope, moment, and shear in bending and amplitudes and phase angles for 
curves of angular deflection and torque in torsion can be computed. The 
phase angles and amplitudes for the deflection and shear curves in bending 
(fig. 2) and the phase angles and amplitudes for the angular displacement 
and torque in torsion (fig. 3) have been computed with reference to a 
unit tip deflection for the weight position = 17 inches. In figure 4 
the amplitudes in deflection and shear in bending from figure 2 are 
compared with the deflection and shear curves due to the fundamental 
uncoupled bending mode of the wing, and in figure 5 the amplitudes in 
angular deflection and torque in torsion from figure 3 are compared with 
the angular deflection and torque curves due to the fundamental uncoupled 
mode in torsion. There is a notable difference in the shape of the 
amplitude curves computed by the present method and those computed from 
the first normal modes . This discrepancy indicates that several modes 
would have to be employed to obtain satisfactory results by the Rayleigh 
type analysis . 


CONCLUDING RIMARKS 


The method discussed in this paper is not limited to a uniform 
cantilever wing with a single weight. By proper attention to the boundary 
conditions the theory can quite easily be extended to apply to a uniform 
wing carrying any number of arbitrarily placed weights, one of which 
might be considered as a fuselage and made to yield the so-called 
symmetric and antisymmetric types of flutter. Furthermore, for conven- 
ience of application, theoretical values of two-dimensional air forces 
have been used. However, since the method does not depend on the 
particular form of air forces involved, any known or available aero- 
dynamic data could be used. In any event, the method is tedious and 
would, therefore, not be recommended over the Rayleigh type analysis 
when it might be known that only the first few normal modes of the 
structure are sufficient to give satisfactory results. 

For wings that are not uniform the differential equations for 
flutter conditions reduce to ordinary differential equations with 
variable coefficients. In this case the solution would, in general, be 
much more difficult to obtain. For general cases there would be no 
advantage in the operational method of solution although an iterative 
process probably might be used to great advantage* 
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In conclusion it is pointed out that the location of such masses 
as engines, landing gears, and fuel tanks might "be used to advantage in 
increasing the flutter speed of a given wing. As shown "by the particular 
problem analyzed herein and by other experiences a definite region exists, 
peculiar to a given wing, in which masses added forward of the elastic 
axis of the wing tend to increase the flutter speed of the wing. 


Langley Aeronautical Laboratory 

National Advisofy Committee for Aeronautics 
Langley Field, Va., November 30, 19^8 


18 


NACA TN No. 1848 


APPENDIX A 


OUTLINE AND EXTENSION OF METHODS OF FLUTTER ANALYSIS 
AS PRESENTED IN REFERENCES 3 AND 4 
1. Derivation of the Differential Equations That Govern 
the Motion of a Wing at Flutter Speed 


Consider a spanwise element of incremental length dx at station x 
of a wing oscillating in "bending and torsion in a free stream of fluid 
(see Bketch) . 



The displacements Y and 0 of an element of the elastic axis are 
functions of x and t. In order that this element remain in dynamic 
equilibrium the external forces and moments on the element must "balance 
the inertia forces and moments. 

The external forces and moments consist of transverse shearing 
forces and torsional moments, which are transmitted from one element of 
the wing to the next, plus the aerodynamic lift force and pitching 
mcment and internal or structural damping. Structural damping is not 
taken into consideration in this discussion, although its inclusion 
would add no computational difficulties . 
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The transverse shearing foroe acting upward at x Is 

N = -EI b ^ (Al) 

cPx 

and that acting downward at (x + dx) is 

N + H ** = " EI b 0 " EI b f“i (A2) 

Similarily the nose-down torsional moment acting at x is 

T = (A3) 

and at (x + dx) the nose- up torsional moment is 

T + s ajt = GJ if + aJ f? dx (M) 

The two-dimensional aerodynamic forces acting on an elemeiit dx of an 
oscillating airfoil have been derived by Theodorsen (reference 8) and 
can be written as a lift force and aerodynamic moment acting about the 
elastic axis of the wing, respectively, a& 


L dx = (o) 2 LyY + oiLy ' || + cd 2 !^© + (4% * ^)dx 


(A5) 


M 


dx = ^o)SlyY + ^ + <uMq 1 -^|^dx 


(a6) 


The inertia force of the element dx can be written 


fm + me! 

V 8t^ 8t / 


and the inertia moment as 


(A7) 
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Diagrams of the forces and moments acting on an element of wing of 
length dx at station x are as follows : 



Imposing the conditions of dynamic equilibrium of the element at x by 
equating inertia forces to external forces and inertia moments to external 
moments gives the two differential equations that govern the motion of 
the wing: 



2 . Boundary Conditions for a Uniform Cantilever Wing Carrying 

an Arbitrarily Placed Weight at Flutter Speed 

The boundary conditions that must be imposed upon equations (A9) 
for a uniform cantilever wing are 

(1) Y(0,t) = 0 

( 2 ) = 0 
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\ 


(3) 0(0, t) = 0 


= 0 


= 0 


0 


that must he Imposed on a vibrating canti- 
lever beam. Condition (l) is the condition that the end at x = 0 is 
supported (either hinged or built in). Conditions (2) and ( 3 ) imply 
that this end is fixed or built in. Conditions 00, (5), and (6) imply, 
respectively, that there is no bending moment, transverse shearing force, 
or torsional moment acting at the tip x = Z. 

If there is an arbitrarily placed weight on the wing, other condi- 
tions must be imposed that will determine the effect of the weight upon 
the motion of the wing. If the weight is considered as concentrated at 
Borne point on the chord line at station x = Z^, it will create discon- 
tinuities in both transverse shear and torsional moment. The magnitude 
of these discontinuities are known functions of the mass of the weight, 
the location of the weight, and the acceleration of the wing. The 
remaining conditions required to complete the boundary- value problem for 
the general motion of the weighted wing are, therefore. 


00 EI t 
(5) 


£ T(x,t i 


x=2 


i& H 


(6) GJ 


Si e(I ' t) 


d=2 


These are the usual conditions 


(7) EI> 


fel Y(x t) 

Y(x,t) 


b z ’Jx.( V 0) 

8x3 

x=(li+0\ 


V 

g 


f? T(x ' tJ * * gi 


( 8 ) 


J z=(Z 1 -0) 


’I” 0(x, t)| 

- -*x=( z 1 +0) 


w. 


w 


82 


2 ^2 


si d 8t 2 


^-Y(x,t) + K 2 ^ 0 (x,t) 


at* 


— * x=2j_ 
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For the purpose of flutter analysis it is assumed that the motions 
in "both bending and torsion are harmonic and that the frequencies in 
tending and torsion are egual. Therefore, only the particular form that 
the solution to the "boundary- value problem has when these conditions 
obtain need be sought. These conditions imply that Y(x,t) and 0(x,t) 
are of the forms 

Y(x,t) = y(x)e tot 

. . f (Aid 

©(x,t) = 0(x)e 

where, on the right-hand side of equations (A10), y and 0 are now 
complex amplitude functions of the span coordinate x from which the 
shape and phase relation of the wing at any fixed time during flutter 
can be obtained. 

If the values of Y and 8 from equations (A10) are substituted 
into both differential equations (a 9) and into the boundary conditions, 
the problem is greatly simplified. The differential equations become 
independent of t and appear as ordinary differential equations with 
constant coefficients . After making the substitution and rearranging 
terms, the equations of motion can be written as 


Elb - ^m + Ly + iLy o£y - ^mep + I# + iLe 0^0 = 0 
GJ ~ + ^me-L + My + iMy + (^1 + M@ + iMg 'J a£e = 0 


(All) 


or more simply as 



ay - 00 = 0 


> 


— — + yy + 50 = 0 
dx 2 


(A12) 
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The boundary conditions also become independent of t and can he written 
as follows: 


( 1 ') 

<< 

O 

II 

O 





(2 1 ) 

y',(o) = 0 





(3') 

e(o) = c 





(V) 

y" 

( Z) =0 





(5') 

y" 

'(I) = 0 





(6 1 ) 

e'(z) = 0 





(7') 

Elfc 

y"'(*i - 

0) - y' 1 'Oi 

+ 0 ) 

II 

"I** 

PP 

C-J 

H 

w 

+ 02 © (Zi) 

(8’) 

GJ 

0*(Zi - 0) 

- e’(U + 0 ) 

II 

“I** 




3 « Solution of Boundary -Value Problems in Ordinary Differential 

Equations by Operational Methods and Application to a Beam 
Carrying an Arbitrarily Placed Weight 

The boundary- value problem given by equations (AI2) and conditions (l T ) 
to (8 1 ) can be solved by straightforward methods of solving ordinary 
differential equations with constant coefficients. The operational 
method, however, is a much easier and shorter approach, particularly in 
view of the discontinuities in shear and torque. 

Briefly, the solution of a boundary- value problem by operatio na l 
methods consists of applying the Laplace transform to the differential 
equations, the initial conditions (root conditions when applied to beam 
problems), and certain forms of other boundary conditions } of solving 
the resulting system for the transform of each dependent variable j and 
then by applying the inversion integral to the results. The remaining 
boundary conditions are then used to set up relations among whatever 
undetermined parameters that might remain. 

In the case of flutter analysis a complete solution to the equations 
is not needed but only the conditions under which an unstable equilibrium 
may exist. The relations that can be set up between the undetermined 
parameters correspond precisely to this condition. In other words these 
relations appear as a system of homogeneous equations and the satisfaction 
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of the condition that this system of equations have a common solution 
other than the trivial solution corresponds to the "border-line condition 
separating the damped and undamped oscillations of the wing. 

The Laplace transform of f(x) is 


= f(s) (A13.) 


,jf(x)j = 


e -sx f ( x ) a* 


where s may be real or complex and x >0. The sufficient conditions 
that this infinite integral exist are that f(x) have no infinite discon- 
tinuities for x ^0 and that f(x) be of exponential order as x — > » . 
(See reference 9*) In other words finite discontinuities such as those 
appearing in the foregoing problem do not invalidate the operational 
approach . 

The Laplace transform of the nth derivative of a continuous 
f motion with continuous derivatives, for which the function and all its 
derivatives are of exponential order, can be obtained directly from 
equation (A13) as 


L-jf n (x)j- = s^s) - s 11 - 1 ^) - s n - 2 f'(0) f^O) (Al4) 


The Laplace transform is linear in the same sense as differentiation 
or integration. That is, if ai and b^ are constants 

ijanf^x) + a n . 1 f n_1 (x) + . . . + a^x) + b^U) + . . . + b o 0(x)| 

= a n L|f n (x)| + a n _ 1 L^f n " 1 (x)y + . . . + a^jf (x)j- 
+ b m L{e m (x)j + . . . + b 0 Lje(x)| (A15) 


Thus the Laplace transform of a linear differential equation with constant 
coefficients is generally a sum of expressions similar to equation (Al4) . 
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In equation (Al4) the quantities f(0), f'(0), . • • , f 11 ( 0 ) 

are the "boundary conditions at the origin of the dependent variable 
(wing root) that corresponds to constants of integration. When these 
quantities are given they are put directly into the transformed equation. 
When the quantities are not given they correspond to what has been 
called undetermined parameters in the preceding paragraphs and must 
later be determined in terms of other boundary conditions. 

Finite discontinuities in a function or any of its derivatives 
are taken into account by proper attention to the limiting values that 
the function or its derivatives have on the two sides of the disconti- 
nuity. In particular. If a function and its first n derivatives are 
of exponential order, if the first (n - 2) derivatives are continuous, 

If the (n - I) 8 '*' derivative has a finite discontinuity at Xq, and if 
the nth derivative Is continuous except for a singular point at Xq, 
(see sketch) 



the Laplace transform of the nth derivative has the form 

Ljf n (x)J = s^f (s) - B n-1 f(0) - . . . 

- sf n " 2 (0) - e- S3C o[f n - 1 (xo + 0) - f n ~Hx o - 0)] (Al6) 

where f(xo + 0) is the value of f(x) as x approaches Xq from the 
right and f(xo - 0) is the value of f (x) as x approaches Xq from 
the left. In other words the terms in the brackets express the magnitude 

of the discontinuity in ^"^(x) at Xq in the (n - l) 8 ^ derivative 

X Q • 
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An examination of the "boundary- value problem, equation (A12) shows 
that the transform will be given by a sum. of expressions precisely of 
the form of e quation (Al6) . 

In order to interpret the transformed function f(s) in terms of 
the original function f (x) , use may be made of the inversion integral 
discussed in text books on operational calculus j or one may refer 
directly to tables of transform- 

As a Bimple example the operational method is applied to a canti- 
lever beam carrying an arbitrarily placed weight and assumed to be 
vibrating in a vacuum in bending only. 

The boundary-value problem for this case can be written 

EI b^f = (A17) 


(a) 

y(o)= y ' (o) =o 

0>) 

y"(D - j'"U) 

(c) 

El-b 

- 0) 


= 0 

- y ,n (2i + 0 ) 



> 

I 


g 


y(2i) J 


(A18) 


where the symbols have the same meaning as in equation (A12) . 

If the root conditions (a) and the boundary condition (c) are used, 
the transformed problem solved for y(s) giveB 


F(s) 


SY 2 Y 3 V vM 

J - a k s 1 * - a 4 a*- a 4 


e~ Bl l 


(A19) 


where, for brevity, Yg 


y M (0), Y 3 = y"’(0), and 


k m r 

a m m 
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The inverse transform of equation (A19) is (see pair nos. 31 
and 32, p. 296, and relation 22, p. 29b , of reference 9) 


. . Y 2 ( \ Y 3 ( \ 

y(x) = — -IcoBh ax - cos ax )+ — sinh ax - Bin ax I 
2or\ J 2aJ , / 


*/* , . 
+ — o — y(*i) 
2a^gEI- b 


sirih a(x - Zqj - sin 


a(x - li)J 


(A20) 


or 


y(x) = cosh ax - cos ax\ + — sinh ax - sin ax) 
2 af\ ) 2 a 3 \ / 


^cosh aZj - cos aZ^ 


p 

W^CD 

Y 2 

20 ^ 

2 a 2 


l 3 / 

+ sinh aZi - sin aZ^ 


2 a. 3 \ 




sinh a(x - Zj) 


sin a(x - Zi) 


(A21) 


where the last "bracket is zero when x - 2i < 0 . 

Imposing "bo undar y conditions (b) gives two homogeneous equations 
in Y 2 and Y3. Each value of a that will cause the determinant of 
the coefficients of Y 2 and Y3 to vanish corresponds to a mode of 
vibration. 

This result has been applied to the wing and weight discussed in 
the text of this paper with the weight located 17 inches from the root. 
The deflection and shear curves due to the first uncoupled modes in 
bending only have been computed and are plotted in figure 6. Corre- 
sponding results have been computed by a 20-station process of iteration 
discussed in reference 12 and plotted in the same figure. 
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4. Representation of the Inverse Transform of the Boundary - 
Yalue Problem, Equation (A12), by a Power Series 
The transform of both y(x) and e(x) of equation (A12) are of the 

form 


f(s) 


Pl(s) P 2 (s) 

^r + ‘irsr 


e -I o e 


(A22) 


where Pi(s) and P 2 (s) are polynomials both of lower degree than q(s) . 
Neither Pi(s) or P 2 (s) have common factors with q(e) where q(s) is 
of the specific form 


q(s) = s^ + as^ + bs 2 + c = (s 2 - Ri)(s 2 - " R 3 ) (A23) 

where the coefficients a, b, and c and the roots squared R^, Rg, 
and P 3 are complex. The inverse function associated with such a 
transf orm gives f (x) in terms of circular and hyperbolic functions 
of XnJ Ri, but with the results in this form the process of solving the 
flutter determinant bee ernes very cumbersoma. 

By making use of the properties of symmetric functions. Go land and 
Luke (reference 4) outlined a simple method of obtaining series expansions 
for the transforms of equations (A12) that does not involve the 
meticulous task of finding the roots of q(s) * The inversions of these 
expansions give y(x) and 0 (x) in the form of convergent series. 

For the development of these series it is first necessary to con- 
sider q(s) as a cubic in s 2 } namely. 


,(„) - g^s 2 - Bi) - - £) 


(A24) 


By making use of the 'binomial theorem^ l/q(s) can he written as 


1 1 3 / % Jl 2 Ji 3 \ 


(A25) 
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Equation (A 25 ) is independent of any interchange of the parameters Rp, R 2 , 
and R 3 and thus satisfies the description of a symmetric function in 
these parameters. (For a discussion of symmetric functions see refer- 
ence 10 or any text on higher algebra or theory of equations.) If the 
indicated multiplication in equation (A 25 ) is carried out, the results 
can he written 





+ 



(A2 6) 


where the general term T n represents the sum of all possible symmetric 
polynomials in Rp, Rg, and R 3 which are of degree n and with all 
coefficients unity. By making use of Newton' s identity relative to 
symmetric polynomials, that is 


T n - “&T n _p ” ^n-2 ~ cT n ~3 


(A27) 


where the value of any T n _ . is to be disregarded when n - J <0, 
every Tn can h e written in teims of the coefficients a, b, and c of 
equation (A23)j for example. 



12 = a 2 - b 

T^ = -a^ + 2 ab - c 


(A28) 


With the aid of equation (A26) and equations (9) and (l0)_of the 
text the inverse transform of equation (A22) or of y(s) and 6(b) can 
therefore be written as a sum of terms of the type given in equations (9) 
and (10) where the T n 's enter as coefficients in the numerator and are 
easily evaluated in terms of the coefficients of a known cubic equation. 
In the application to flutter analysis only the first few T n 's are 
usually necessary because the resulting series is generally found to be 
highly convergent. 
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APPENDIX B 

DERIVATION OF THE FLUTTER DETERMINANT AND SAMPLE CALCULATIONS 

Introduction 


In this section the flutter determinant is formally derived and 
the method described in the text for solving the determinant is 
illustrated with sample calculations for a specific example . Also final 
expressions for the deflection curves are given from which amplitude and 
phase angle curves of deflection, shear, and torque are calculated for 
a specific case. The calculated amplitudes are compared with corre- 
sponding curveB computed from the fundamental uncoupled modes In bending 
and torsion. 


Derivation of the Flutter Determinant 

In equations (ll) and (12 ) of the text It Is first necessary to 
evaluate the expressions 


- 0) - y’"Ul + of 


and 


e 


'(ii - o) - 



in terms of Y2, Yo, and 0^. Since terms involving (x - Ip) drop out 
of both equation (ll) and equation (12) for x = Zp, the values of y(l^) 
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nnfl 0(Zi) .'.an "be obtained directly from these equations. The values 
of y(li) and 0(li) substituted into conditions (c) and (d) of the text 
give the desired relations; namely. 


y ,M (Zl - 0) - y'^dx + 0) 




brCii) 


+ 62© (ll) 


EI bS 
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and 

0'Ul - 0) - e'(2i+ 0) = ^y(Zi) + Kg 2 0(Zi)] 



Substituting equations (Bl) and (B2) into equations (ll) and (22) 

gives 

j(x) = hi(x)Y 2 + h 2 (x)Y 3 + h 3 (x)0 1 (B3) 

(B4) 


9(x) = glU)Y2 + 82( x ) y 3 + 83 W 0 ! 
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where 


. . \ T n x? n+2 'T -52 ' T n ^ n+ ^ 

h x (x) =2^ ( 2n + 2 )i + ) (2n + 4) I + El^g 


2n+4 


n=<- 


h=0 


T n li 2n+e 
(2n + 2) i 


IT. 7 

(6 -ejr) 2_ -^-rtTT 

- n=0 


T n(? - 0 2n+5 , r 1 ! 1 - 0 2n+3 

4 * + 4 o < 2n+3 > ! - 


fen +3 )! J 


PVfyj£ 


GJg l: 


" T n Z! 2n+4 


(e 2 B - *s 2 ) 21 + It) I 


n=0 


m 7 2n+2 — | 
T n z l I 


V~ T n t l \ T n( z " 4 lJ 

^^bTSTTsJ tJ 


00 T n (x - Z!) 2n+5 


T n ^ n+3 . «^I 7 , 


\ °° T n X? n+ 5 — T n X-“ -’ 

h 2 u> = 6 2_ (2n + j, + 2_( 2 „ + 3 )! + 
n =0 n =0 


EI hS 


T n Z 1 2n+5 
(5- e27)}_ (gnt ' sP. 

- n=0 


, \^ T nA 2nt3 ~ 

4^TS— srr 


u. 


X - 


(2n + 5) I 


T n( ]r ~ hf** 3 ' 
fen + 3) 1 


PV^ 

GJg 


» T 7 2n+5 
(625 - Kj2 2 7) 2 _ ( gn + 5 ) ■ 


n=0 


3 ^ 
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h 3 ( x ) 


81 (x) 


82 ( x ) = 


„r v 211 ^ 

7 ^TTTVT Ti!T r * 




00 


+ 


e 2 



2n+l“ 


(2n + 1) i 




f* - ■*> 



W a ? 

v 

ei T )8 


v — — — m 7 2 n +4 — m 1 2n+2 

\ Vi ^ 6 i 

Z— (2n + i|) I Z (2n + 2) ! 
n=0 n=0 


f^Tjx - 

fen + 5) ! 


n=0 


TnX 2 n +5 


_ \ ' 

7 Z _ fen + 5 ) I + GJi" 

n=0 


V" Vl 2n+3 ' 

- *2 Z _ tstts 3 '! 

n=0 


a 


7 2 \ T n H 2n+ 5 

(9BS ' K2r) II (gn\ ? ); 

— n=0 

T n(* - >l) 2n+5 


n=0 


2n + 5)1 


f- T n ( x - 2 1 )2n+r 
Z _"" 72n”iyi — _ 


n =0 


W ^7 


EI bS L_ 


(6 


w 

527 ) 


11=0 


Vi 2n+ ^ 

(2n + 5) I 


rvi 2n+3 v T n ( 

2__ tst in z_~~ 

n=0 n=0 


TJ x - l ! 


.2n + 5, 


) 2 n +5 

n — 
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and 


g 3 (x) 


V 2 - T nX 2 n+1 v 2 - T n x 2n+ 5 

^TSTTTTT - 5): 




GJg H 


(e2P - Kj2 2 a) 



TnU 2n+? 
(2n + 5)1 


*2 


, 2 > 




T Zi 2n+ ^ 
-‘•n'-l 

' f : 

TnU - H ) 2n+5 C 

• T n U - 

(2n + i)\] 

- « 

(2n + 5 ) • 

j (2n + 1) ! 


EI t g 


00 

(P - e 2 a) 'y~ 


?n?l 2n+ 5 
(2n + 5) I 


n=0 


+ 


N^~ T n 2-L 2n+1 

e 2 Z_ (2n + l) I 
n=0 



By imposing conditions (t>) of the text 

y"(l) - y MI U) = e'O) =0 

upon equations (B3) and (b 4), three equations are obtained (written in the 
text as equation (l3)) : 


A ± Y 2 + BiY 3 + Ci©! = 0 
where i = 1, 2, and 3 and 


A x = ^"(z) 

Ag = hj' ' 1 (l) 

a 3 = gi'(i) 


B x = hg M (Z) 

Bg = hg'^Z) 
B 3 = efe '(2) 


C x = h 3 "(Z) 
Cg = h 3 "'(Z) 

c 3 = g 3 1 ( l) 
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Imposing the condition that the equations (13) have a solution 
other than the trivial solution Y 2 = Y 3 = 0 ± = 0 results in the flutter 
determinant 

A i B i c i 

A = Ag Bg C 2 =0 (B5) 

a 3 B 3 C 3 

Sample Calculation of Flutter Speed 
and Deflection Curves 

A method of solving the flutter determinant given in the text is 
illustrated here by the solution of the determinant for the wing-weight 
combination discussed in the text when the spanwise location of the 

weight is 17 inches from the root. The values of — = 7= that are 
^ to lr 

chosen are in the neighborhood of the experimental value and have 

available tabulated values of Theodorsen's function C(k) = F + iG. 

Table I shows the actual computations required to evaluate the 
coefficients Aj^ B^, and C^ for ~ = 7*1429 (k = 0.l4) and two 
values of ~ » f (f = 25 cps and f = 28 cps) . From columns (32) , (38) , 
and (5+1) the determinant for f = 25 cps is 

(14.9200 - 2.85741) (12.8320 - 2.03151) -(7*3286 - 0.600211) 

A = (11.8000 - 3*66951) (10.2970 - 2.85661) -(5.4711 - 0.93233i) 

(0.17030 - 0. 661341) -(0.09077 + 0. 593^11) -(0.41138 - 0.288641) 

or 

A = 1.0326 - 0.69^1 
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Similarly, for f = 28 cps. 


A = 


(18.6380 - 3.81151) 
(15-5930 - 5-09351) 
-(0.04177 + 0.870981) 


(15.0860 - 
(13.OO80 - 
-(0.23526 + 


2.63991) 

3-79461) 

0.759**8i) 


-(9.1238 - 

-(7.1158 - 
-(0.51403 


0.854331) 

1.39881) 

- 0.370171) 


or 


A = -0.4029 - 0.03121 


The determinant was evaluated in this manner for the same value 
of v/bcu and several other values of f. The process was then repeated 

i" 01 * = 8.25 and several values of f and for = 5-00 and several 

values of f . The real and Imaginary parts of "the evaluated dete rmin ant 
for each value of v/bcu and the corresponding values of f are separately 
plotted In figure 7 • The ordinates of the Intersections of the different 
pairs of curves of real and Imaginary parts were scaled in figure 7 and 
plotted as A 3 against both v/bcu and f in figure 8. The zero 
ordinates of these curves give the value of v/bcu = 6.93J and the 

values of f(f = 28.04 cps) for which the determinant vanishes. From 
these values the flutter speed is readily calculated to be 


v = (to) (6.93) = (2nbf) (6.93) = fe«H 2 g- Ql 0( 6 -93)' = 407 fps 


As pointed out in appendix A the deflection curves at any specified 
time are given by equations (AlO) 


Y(x,t) = y(x)e : * m ^ = y(x)(cos cut + i sin cut) 
0(x,t) = ©(xje^^ = 0(x)(cob cut + i sin cut)^ 


where final forms of y(x) and 0(x) are given by equations (B3) and (B4) 
and where, at least, the relative values of the undetermined coeffi- 
cients Y2, Y3, and 0p in equations (B3) and (B4) must be known. If 
the set of values of v/bcu and cu that satisfy the flutter determinant 
is used to determine the coefficients A^, Bj, and in equations (13)> 

there is obtained a system of three homogeneous equations in the three 
unknowns Y2, Y3, and 6 -\_ that have solutions other than the trivial 


solutions Y^ = Y^ = 0^ = 0. If these equations are each divided 
through by any one of the unknowns, say Yg , there is obtained a 
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consistent system of three equations in the two ratios Y^/Y 2 and 0^/Y 2 • 

Any two of the three equations can therefore he solved for these ratios. 
Consequently, equations (B3) and (B4) can he written with one undetermined 
parameter that appears as a factor in each equation. Furthermore , since 
the coefficients B^, and are complex numbers the ratios Y i/Y 2 

and 0i/Y2 are complex numbers and equations (B3) and (B4) contain 
complex coefficients. The real and imaginary parts of these equations 
can he separated and the equations written 

(B 6) 


y(x) = Y 2 [ji(x) + iy 2 (x)j 
0(x) = Y 2 [e 2 (x) + i0 3 (x)] 

If these relations are substituted into equations (AlO) , 


Y(x,t) = Y 2 |yj(x)cos cut - y 2 (x)sin cot + iJ^y 2 (x)cos cot + y^(x)sin cot||] 


1 


(BT) 


e(x,t) = Y 2 |0 2 (x)cos cot - 0 3 (x)sin cot + i|j> 3 (x)cos cot + 0 2 (x)sin cot j 


or 


Y(x, t) =Y 2 


2 

x)| + 


[y 2 (x)f [cos (cot + cpj.) + i sin(u>t + <Pl)] 


e(x,t) = y 2 \ [ 0 2 ( x )] 2 + (0 3 (x)] 2 Jcob (cot + (f^) + i sin (cot + q^>)J 


(B8) 


where 


and 


CPi = tan 


y 2 (*) 

y x (x) 


. -l ® 3 ^ x ^ 
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and where <p-]_ - cpg represents the difference in phase angle between 
bending motion and torsion motions at x. 

The real parts of equations (B8) are interpreted to mean the motions 
in bending and torsion taken in a positive sense* The imaginary parts 
can then be interpreted as representing these same motions with a phase 
shift of 7t /2 radians . 
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Ratio of flutter speed of wing with weight to flutter speed of wing without weight, 
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Distance along span, In. 

Figure 1.— Comparison of calculated and experimental flutter speeds for 

a particular wing-weight system. 
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Dietanoe along span, In. 


Figure 2. Plot, of amplitude and phase angle of displacement and shear 
curve in bending at flutter for = 17 inches (amplitude and shear 
referred to unit amplitude at tip in bending) . 


Phase angle, deg 
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Distanoe along span, In. 


Figure 3*— Plot of amplitude and phase angle of torsional displacement 
and torque for 2q = 17 incheB at flutter (amplitude and torque 
referred to unit amplitude at tip in bending) . 
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Figure 4.— Plots of amplitudes in "bending displacement and torque and 
the corresponding curves computed for the first uncoupled normal mode 
in tending for l± = 17 inches (amplitude and shear referred to unit 

amplitude at the tip in tending) . 
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Figure 5*— Plots of amplitudes in torsional displacement and torque and 
the corresponding curves computed for the first uncoupled normal mode 
in torsion for Z-]_ = 17 inches (amplitude and torque referred to 
unit amplitude at tip in tending) . 
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Figure 6.— Plots of deflection and shear curves computed from the first uncoupled modes in "bending 
by the differential equation method and by the 20— station iteration process of reference 12 
(referred to unit tip deflection) . 







